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(N \ Abstract 



. The ARGO-YBJ experiment is currently under construction at the Yang- 

O \ bajing Cosmic Ray Laboratory (4300 m a.s.l.). The detector will cover 74 x 78 m 2 

with a single layer of Resistive Plate Counters (RPCs), surrounded by a partially 

^ instrumented guard ring. Signals from each RPC are picked-up with 80 read out 

^ | strips 6 cm wide and 62 cm long. These strips allow one to count the particle 

number of small size air showers. In this paper we discuss the digital response 

of the detector for showers with core located in a small fiducial area inside the 

carpet. The results enable us to assess the sensitivity of the strip size spectrum 

measurement to discriminate between different models of the Primary Cosmic 

"35 ', Ray composition in the energy range 10 -r- 500 TeV. 

d • 

1. Introduction 

X 

H ' The energy spectrum of cosmic rays is well described by a power law over 

several decades of energy, before and after the so called knee region, 10 eV -j- 
10 16 eV, where the slope changes. Despite of many conjectures and attempts, the 
origin of this steepening, observed in air shower data, is still obscure. Comparing 
existing data makes evident a substantial disagreement between the primary cos- 
mic ray composition models provided by different experiments. As an example, 
the proton spectrum measured by TIBET AS 1 -/ [1] changes its slope at energy 
around 100 TeV, KASCADE [2] data suggest a steepening at about 2 PeV 
while data collected by the balloon-born experiment RUN JOB [3] dont exhibit 
any spectral break up to 500 TeV . In this respect, ARGO-YBJ offers an unique 
possibility of measuring small size air showers down to a few TeV primary energy 
by exploiting the digital read-out of the detector. In this paper we discuss the 
sensitivity of the ARGO-YBJ detector to discriminate between different models 
of the Primary Cosmic Ray composition in the energy range 10 -r- 500 TeV. 
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Fig. 1. Schematic view of the 
ARGO-YBJ detector . 
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Fig. 2. Average strip size and pad 
size compared to the total size and 
truncated size . 



2. The ARGO-YBJ detector and the digital read-out 

In the ARGO-YBJ experiment (Fig. 1 and [9]) the signals from each RPC 
are picked-up with 80 read out strips 6.7 cm wide and 62 cm long for a total of 
124800 in the central carpet, with an average density of ~ 22 strips /m 2 . The 
FAST-OR of 8 strips defines a logic unit called PAD (10 PADs for each RPC) this 
signal being used for timing and trigger purpose. The detector is clustered in units 
of 12 chambers (CLUSTER) with modular read-out (see Fig. 1). A simulation has 
been carried out by means of the CORSIKA/QGSjet code [7] in order to study 
the dependence on the energy of the number of fired strips for quasi-vertical 
(< 15°) showers with core in a fiducial area (Af) of ~ 260 m 2 at the center of the 
carpet (2x3 CLUSTER). An average strip efficiency of 95% and an average strip 
multiplicity m = 1.2 have been taken into account. The average strip size (N s ) 
and pad size are compared in Fig. 2 to the total size and to the size sampled by 
the central carpet (truncated size) for proton induced showers. The Fig. 2 shows 
clearly that the digital response of the detector can be used to study the primary 
spectrum up to energies of a few hundreds TeV. 

3. The strip size spectrum 

To check the sensitivity of a digital measurement to the primary spectrum, 
the strip size spectrum has been obtained for four models of the primary cosmic 
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ray composition by calculating the differential ffux 

dN Z" 00 

- r = Jaii{y s )=Y. e A{ys)- / J A (E)-P A (y s ,E)-dE (1) 

where y s = Log(N s ), Ja{E) is the differential intensity for each nuclear species, 
P A (y s ,E) is the probability that a nucleus of mass number A and energy E 
produces a strip size y s . The following models for the primary composition have 
been considered: J AC EE [4]; RUN JOB [3]; TIBET ASj [1] and the model 
proposed by J. R. Horandel [6] to fit world data from EAS experiments. To get 
Pa{Vs,E) we have envisaged the following experimental conditions: 1) a trigger 
requiring at last 16 PADs out of 120 in any one of the 4x5 central CLUSTER 
(~ 1000 m?) including the fiducial area Af, 2) shower core and shower direction 
reconstructed following the algorithms developed in [5] and [8] respectively. The 
probalilities for shower with core recostructed in Af and zenith angle 9 < 15° may 
be described by a gaussian form P A (y s , E) = ^ 1 ■ exp[-\ ■ [ ^"ffiff ] 2 ]- 
Showers which dont satisfy the trigger condition or are reconstructed with core 
outside the fiducial area are not considered. The efficiency £a(Vs) , see Fig. 3, 
represents the fraction of triggered and reconstructed events with core in Af. In 
Fig. 4 is plotted the function F all (y s ) = J a u(y s ) ■ 10 1,5 ' ys T- Ay s with Ay s = 0.1 and 
T = 10 s m 2 ssr. This function provides in a simple way the slope of the strip size 
spectrum {J aU {N s ) oc N~ a °) in the interval y n ^y S2 , a s = 2 .5- L ° 9[Fal f g s ^ auM 
and the counting rate number of collected events, C s (y s ), integrated in Ay s = 0.1 
for an exposure T = 10 8 m 2 ssr, corresponding to about one month of data taking 
for showers with core selected in Af and zenith angle 9 < 15°, C s (y s ) = f^'i^ ■ 

4. Discussion and conclusions 

From the results of Fig. 4 we can draw the following conclusions: 

1. the strip size spectra pile-up above y s ~ 4.2 (< E p >~ 100 TeV, 
< Ep e >~ 320 TeV), due to the saturation of the digital read-out, and fall down 
below y s ~ 3 (< E p >~ 10 TeV, < E Fe >~ 40 TeV) owing to a substantial 
decrease of the trigger efficiency (e(y s = 3) ~ 75%, see Fig. 4); 

2. in the range y s = 3 -j- 4 the JACEE data and the Horandel model 
predict the same strip size spectrum with a s ~ 2.35. This is mainly due to the 
fact that in these models the spectrum of light component is quite similar; 

3. the biggest difference is found between the strip size spectra obtained us- 
ing the TIBET AS'y and RUN JOB models. The spectral index is ~ 2.27 and ~ 
2.37, respectively, and the counting rate expected according to the TIBET ASj 
spectrum is higher than the one predicted by the RU N JOB data, running from 
~ 30% at y s ~ 3 up to ~ 55% at y s ~ 4; 

4. the number of events in each size bin is enough to make negligible 
the statistical uncertainty. On the contrary, any systematic error 5N S /N S in 
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Fig. 3. Fraction of triggered and re- 
costructed events with core in the fidu- 
cial area for proton and iron induced 
air shower. 
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Fig. 4. Size strip spectrum for four 
models of primary composition (see 
text). 



reconstructing the strip size spectrum determines a shift AF a u/F a u = (5N S /N S + 
l) as_1 — 1. Thus, a control of the detector performance at a level better than 10% 
is required in order to reduce any systematic effect below 15%; 

5. the spectra expected from using a different hadronic model have not 
been calculated. Due to low energy range involved, we expect differences lower 
than 10%. 

The possibility of discriminating between different models of the primary 
cosmic ray composition by using the digital read-out of ARGO-YBJ has been 
investigated. It is important to notice that this technique should allow to scan 
the energy range from 10 to a few hundred TeV where direct and indirect mea- 
surements partially overlap. 
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